V 

AO-A102  626 

UNCLASSIFIED 


NAVAL  AIR  DEVELOPMENT  CENTER  WARMINSTER  PA  AIRCRAFT  —ETC  F/G  1/3 
THE  CONTROL  ANTICIPATION  PARAMETER  FOR  AUGMENTED  AIRCRAFT. (U) 

MAT  81  D  E  BISCHOFF 

NADC-81186-60  _  mi 


ADA102626 


THE  CONTROL  ANTICIPATION  PARAMETER 
FOR  AUGMENTED  AIRCRAFT 


D.  E.  BISCHOFF 

Aircraft  and  Crew  Systems  Technology  Directorate 


NAVAL  AIR  DEVELOPMENT  CENTER 
Warminster,  Pennsylvania  18974 


Approved  for  Public  Release;  Distribution  Unlimited 


Prepared  for 

NAVAL  AIR  SYSTEMS  COMMAND 
Department  of  the  Navy 
Washington,  D.  C.  10361 

81  8  10  049 


NOTICES 


REPORT  NUMBERING  SYSTEM  -  The  numbering  of  technical  project  reports  issued  by  the 
Naval  Air  Development  Center  is  arranged  for  specific  identification  purposes.  Each 
number  consists  of  the  Center  acronym,  the  calendar  year  in  which  the  number  was 
assigned,  the  sequence  number  of  the  report  within  the  specific  calendar  year,  and 
the  official  2-digit  correspondence  code  of  the  Command  Office  or  the  Functional 
Directorate  responsible  for  the  report.  For  example:  Report  No.  NADC-78015-20 
indicates  the  fifteenth  Center  report  for  the  year  1978,  and  prepared  by  the  Systems 
Directorate.  The  numerical  codes  are  as  follows: 


CODE  OFFICE  OR  DIRECTORATE 

00  Commander,  Naval  Air  Development  Center 

01  Technical  Director,  Naval  Air  Development  Center 

02  Comptroller 

10  Directorate  Command  Projects 

20  Systems  Directorate 

30  Sensors  &  Avionics  Technology  Directorate 

40  Communication  5  Navigation  Technology  Directorate 

SO  Software  Computer  Directorate 

60  Aircraft  6  Crew  Systems  Technology  Directorate 

70  Planning  Assessment  Resources 

80  Engineering  Support  Group 


PRODUCT  ENDORSEMENT  -  The  discussion  or  instructions  concerning  commercial  products 
herein  do  not  constitute  an  endorsement  by  the  Government  nor  do  they  convey  or 
imply  the  license  or  right  to  use  such  products. 


% 


m 


■SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whin  D#l.  Cnl»,md) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


H  REPORT  NUMBER 

NADC-81186-60 


2.  GOVT  ACCESSION  NO. 


3  RECIPIENT’S  CATALOG  NUMBER 


/dal 


t 


4  a.  -^TrTLE'i'APd  Subtitii)  ' 

I’  The  Control  Anticipation  Parameter  for 
‘  Augmented  Aircraft . 


S.  TYPE  OF  REPORT  t  PERIOD  COVERED 

Final  1  Jan-»  15  May  81 J 


«.  PERFORMING  ORG.- REPORT- HUMMM • 


7.  AUTHORfiJ  _ 

1  D.  E./Bischoff 


a.  CONTRACT  OR  5BAMT  HUM  BERf  t) 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Aircraft  &  Crew  Systems  Technology  Directorate 
Naval  Air  Development  Center 

Warminster,  PA  18974  _ _ 


10.  program  element,  project,  task 

AREA  S  PORK  UNIT  NUMBERS 

AIRTASK  A034-0000/001B/ 
OF41-400-000 


II.  controlling  office  name  and  adoress 
Naval  Air  Systems  Command 
Department  of  the  Navy 
Washington,  DC  20361 _ 


a 


12,  REPORT  OATS 

15  May  1981 


IS.  NUMBER  d?*"^AGES 

37 


14.  MONITORING  AGENCY  name  »  ADDRESS^!/  dltlbrent  I fom  ConrrolllnJ  Olllcm) 


:  W 


15.  SECURITY  CLASS,  (ol  thl a  report) 


UNCLASSIFIED 


1S«.  OECLASSl  FI  CATION/  DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  this  Report) 


17.  DISTRIBUTION  STATEMENT  (of  the  ebetrect  entered  in  Block  30,  it  different  from  Report) 


Approved  for  Public  Release;  Distribution  Unlimited 


18.  SUPPLEMENTARY  NOTES 


19-  KEY  WORDS  ('Continue  on  reveree  tide  if  neceeemry  end  Identify  by  block  number) 

Aircraft 

Flying  Qualities 
Control  Systems 
Equivalent  Systems 


ABSTRACT  (Continue  on  rover ee  elde  if  neceeeery  end  Identity  by  block  number) 

A  longitudinal  control  anticipation  parameter  is  developed  which  correlates 
high  order  augmented  aircraft  dynamics  with  both  lower  order  equivalent  svstem 
descriptions  and  pilot  opinion  ratings.  Three  extensive  in-flight  research 
programs  are  analyzed  via  the  developed  parameter  in  an  effort  to  define 
regions  of  pilot  acceptability.  Additional  effort  is  recommended  to  verify  the| 
minimum  level  1  boundaries  and  to  determine  the  relationship  between  the  con¬ 
trol  anticipation  parameter  and  alternate  proposed  criteria.  A 


DO 


FORM 
I  JAN  73 


1473  EDITION  OF  I  NOV  ««  IS  OBSOLETE 

S/N  01 02-LF-01 4-4601 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  RAOE  (Whm  Dud  toil* bd) 


NADC-81186-60 


TABLE  OF  CONTENTS 


LIST  OF  FIGURES  . 

LIST  OF  TABLES . 

LIST  OF  SYMBOLS  . 

INTRODUCTION . 

CONTROL  ANTICIPATION  PARAMETER  DEVELOPMENT 

CAP'  BOUNDARY  DETERMINATION  . 

EQUIVALENT  SYSTEM  REPRESENTATION  OF  CAP' . 

FLIGHT  TEST  CORRELATION  OF  CAP'  . 

SUMMARY  . 

REFERENCES . 


2 

2 

3 

4 
1 

17 

17 

24 

36 

37 


1 


NADC-81186-60 


LIST  OF  FIGURES 


Figure  No.  Title  Page 

1  MIL-F-8785  Short  Period  Frequency  Requirements .  5 

2  F-14  Equivalent  Short  Period  Frequency  Parameters  ....  6 

3  Attenuation  Effect  of  Feel  System  Dynamics  on 

Pitch  Response . 10 

4  Attenuation  Effects  of  Aircraft  Dynamics  on  Maximum 

Pitch  Acceleration.  . . 11 

5  Relationship  between  CAP  and  Pilot  Opinion . .  18 

6  Relationship  between  CAP'  and  Pilot  Opinion  .  19 

7  Specification  Requirements  Based  on  CAP  and  CAP* .  20 

8  Equivalent  System  Comparison  -  A-6  Airplane  .  22 

9  Control  Anticipation  Parameter  Comparison  -  NEAL-SMITH 

Data . 29 

10  Control  Anticipation  Parameter  Comparison  -  LAHOS  Data.  .  30 

11  Time  Delay  versus  CAP'  -  Neal-Smith  Data .  ....  32 

12  Time  Delay  versus  CAP'  -  LAHOS  Data . .  33 

13  Bandwidth  Criterion  -  Neal-Smith  Data  . . 34 

14  Bandwidth  Criterion  -  LAHOS  Data . 33 


LIST  OF  TABLES 


Table  No.  Title  Page 

I  -  NEAL-SMITH  Data  Configuration  Summary  .  .  25 

II  LAHOS  Data  Configuration  Summary.  ...  .  26 

III  Parameter  Tabulation  -  NEAL-SMITH  Data . 27 

IV  Parameter  Tabulation  -  LAHOS  Data  .....  .  .  .  28 


2 


NADC-81186-60 


LIST  OF  SYMBOLS 


CAP 

CAP' 

F 

Ke 

La 

n/a 

nz 


s 

td 

T92 

V 

^6e>  ^ 


a 

6c 

6e 

9 

9,q 

9 


Control  Anticipation  Parameter  -  rad/sec^/g 
Attenuated  Control  Anticipation  Parameter  -  rad/sec^/g 
Pilot  force  input  -  lb 

Gain  of  pitch  rate  to  control  input  transfer  function 
Inverse  numerator  time  constant.  Dimensional  lift  curve  slope 
Acceleration  sensitivity  parameter  -  g/rad 
Normal  acceleration  -  g's 

Pitching  moment  stability  derivative  due  to  control  input,  vertical 

velocity  and  pitch  rate,  respectively 

Laplace  frequency  domain  variable 

Equivalent  system  time  delay  -  sec 

Pitch  rate  numerator  time  constant  -  sec~l 

Freestream  velocity  -  ft/sec 

Normal  force  stability  derivative  due  to  control  input  and  vertical 

velocity,  respectively 

Angle  of  attack  -  rad 

Pilot  control  input  -  rad 

Elevator  (stabilizer)  position  -  rad 

Pitch  angle  -  rad 

Pitch  rate  -  rad/sec 

Pitch  acceleration  -  rad/sec2 

Undamped  natural  frequency  -  rad/sec 

Damping  ratio 


Subscripts 


e 

FS 

sp 

S3 

nd 


-  Equivalent  system 

-  Feel  system 

-  Short  period 

-  Steady  state 

-  Nondimensional 


NADC-81186-60 


THE  CONTROL  ANTICIPATION  PARAMETER  FOR 
AUGMENTED  AIRCRAFT 

INTRODUCTION 

The  military  flying  qualities  specif ication,  MIL-F-8785B  (reference  (a)) 
specifies  airplane  short  period  frequency  (wnsp)  requirements  as  a  function  of 
acceleration  sensitivity  (n/u)  as  shown  in  Figure  1.  The  boundaries  for  mini¬ 
mum  and  maximum  frequency  requirements  have  been  established  from  piloted 
flight  tests  as  lines  of  constant  wnSp/n/a.  Using  the  assumption  that  the 
responses  of  interest  to  a  pilot  during  a  pullup  are  the  initial  pitch  acceler¬ 
ation  and  steady  state  normal  acceleration,  Bihrle  (reference  (b))  defined  a 
control  anticipation  parameter  (CAP)  as  the  ratio  of  these  two  parameters.  He 
further  showed  that  in  the  short  period  approximation,  CAP  is  equal  to  mns2/n/ci. 

Both  the  8785B  requirements  and  Bihrle’ s  analysis  were  developed  using 
unaugmented  aircraft  for  which  the  control  system's  dynamic  effects  were  assumed 
to  be  negligible.  If  higher  order  dynamics  are  included  in  the  system  descrip¬ 
tion,  this  simple  relationship  may  no  longer  hold.  For  example,  the  implemen¬ 
tation  of  a  washed  out  pitch  rate  feedback  does  not  alter  either  the  initial 
pitch  acceleration  or  the  study  state  normal  acceleration,  thereby  not  affecting 
CAP.  However,  the  short  period  frequency  may  be  significantly  changed  causing 
u)ns2/n/a  to  differ  from  CAP.  If  feel  system  or  actuator  dynamics  are  included 
in  the  system  model,  the  initial  pitch  acceleration  is  identically  zero, 
building  to  a  maximum  some  time  after  the  input  is  applied.  Difranco  studied 
this  problem  and  defined  on  attenuated  CAP,  CAP',  which  took  into  account  the 
effects  of  the  feel  system's  dynamics,  reference  (c). 

The  latest  revision  to  the  military  flying  qualities  specification, 
MIL-F-8785C  (reference  (d))  imposes  the  short  period  requirements  of  -8785B 
on  aircraft  possessing  higher  order  control  system  dynamics  by  allowing  the 
determination  of  equivalent  lower  order  system  characteristics  which  approximate 
the  higher  order  system.  While  it  may  be  argued  that  it  is  possible  to  determine 
equivalent  lower  order  systems  having  similar  time  responses  to  those  of  the 
actual  higher  order  systems,  correlation  of  the  resulting  parameters  with  MIL- 
SPEC  requirements  has  produced  some  perplexing  results.  For  example,  aircraft 
possessing  complex  control  systems,  whose  dominant  roots  appear  to  lie  in  the 
acceptable  regions  of  Figure  1,  yield  equivalent  systems  whose  versus  n/a 
relationship  crosses  the  minimum  frequency  boundaries  into  unacceptable  regions 
(Figure  2).  However,  the  time  histories  for  such  responses  show  little  differ¬ 
ence  from  those  of  the  higher  order  systems  from  which  they  were  developed. 
Further,  analysis  of  the  control  anticipation  parameter  (by  measuring  the 
maximum  slope  of  pitch  rate  and  steady  state  nz)  for  each  of  these  responses, 
does  not  correlate  either.  Therefore,  it  is  difficult  to  interpret  what  a 
pilot's  opinion  of  such  responses  would  be  from  analysis  of  these  model  para¬ 
meters. 

The  purpose  of  this  memorandum  is  to  provide  a  definition  of  the  control 
anticipation  parameter  for  higher  order  systems  which  will  (1)  correlate  with 
the  frequency/acceleration  sensitivity  relationship  of  the  higher  order  system 
and  (2)  be  consistent  with  the  modal  parameters  obtained  from  the  lower  order 
equivalent.  This  will  be  accomplished  by  briefly  reviewing  Bihrle's  and 
DiFranco's  developments  and  then  extending  them  to  general  higher  order  cases. 
Examples  from  current  Navy  fleet  aircraft  and  contractor  flight  research  pro¬ 
grams  will  be  utilized  to  illustrate  the  analysis. 
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CONTROL  ANTICIPATION  PARAMETER  DEVELOPMENT 


Under  the  assumption  of  constant  speed  equations  of  motion,  the  approximate 
transfer  functions  relating  the  short  period  pitch  response  of  an  aircraft  to 
an  elevator  input  can  be  expressed  as  - 


e<s)  _  Ke  U  -v  ) 

sC^2,  ■+•  ^n^p) 


QCs)  _  s  Q(s"> 
ScCs'i  “  S«C  s) 


Qis)  _  szG(s) 
SeCS)  -  Se.Cs) 


The  initial  pitch  acceleration  response  to  a  unit  step  elevator  input  is 
obtained  from 
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The  steady  state  normal  acceleration  may  be  determined  from  the  relationship 
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Bihrle's  control  anticipation  parameter,  defined  as  the  ratio  of  initial  pitch 
acceleration  to  steady  state  normal  acceleration,  may  be  expressed  as: 


It  can  also  be  shown  from  the  equations  relating  n2  and  a  to  elevator  inputs 
that : 

n  _  <\r\  _  V  / 

01  a*  “  ^  y 


which  for 

reduces  to 


DiFranco  expanded  this  analysis  to  include  feel  system  dynamics.  In  this  case, 
the  pitch  response  to  pilot  force  inputs  is  defined  by: 


MSe  - 
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The  initial  pitch  acceleration  response  to  a  unit  step  force  input  is  now 
found  to  be  zero: 


ew  -  *  -  =  7717 
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o O 


s  s^ecs) 

FCs^> 


r  O 


The  maximum  pitch  acceleration  occurs  at  some  time  (greater  than  t  =  0+)  after 
the  input  is  applied  and  is  attenuated  from  that  obtained  without  feel  system 
dynamics  (Figure  3).  An  attenuation  factor,  9n(j,  relating  the  maximum  pitch 
acceleration,  including  feel  system  dynamics,  to  the  initial  pitch  acceleration, 
excluding  feel  system  dynamics,  in  response  to  a  step  stick,  force  input,  can 
be  expressed  as: 


Q wax  PS 


_ (s  •»  '/TeJ _ 

+  ‘*,r^p)(S2‘+2/f.sa3pSS 


n) 


where  9n<j  is  now  a  nondimensional  pitch  acceleration.  The  value  of  encj  can 
be  determined  by  converting  equation  (9)  to  a  time  function  and  finding  its 
maximum  value.  For  any  particular  aircraft  flight  condition,  as  -osp  is 
increased,  the  value  of  9n(j  will  be  reduced,  i.e.,  the  attenuation  will 
become  greater  (Figure  A) . 
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FIGURE  3.  Attenuation  Effect  of  Feel  System  Dynamics  on  Pitch  Respo 
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FIGURE  4.  Attenuation  Effects  of  Aircraft  Dynamics  on 

Maximum  Pitch  Acceleration  (Reference  (c)  figure) 
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The  steady  state  normal  acceleration,  following  a  step  stick  force  input, 
is  now  determined  as: 


s$ 


a  _  V_  s9(Q  j_ 
1  **s  "  ^  s 


V  tOps  1 

^  cJPS 


(VO) 


which  is  unchanged  from  the  case  in  which  feel  system  dynamics  were  excluded. 
Substituting  equations  (9)  and  (10)  into  the  expression  for  CAP  yields 


CAP 


PS- 

Qnd  ^ss 


Z 

n/ot 


DiFranco  further  defined  CAP'  to  include  the  attenuating  effects  of  feel 
system  dynamics  as 


cap' 


n, 

*  IS 


z 


Wot 


Examination  of  this  last  equation  provides  insight  into  the  problems  of  mapping 
higher  order  system  short  period  characteristics  onto  specification  requirements 
and  attempting  to  correlate  them  with  the  pitch  and  normal  acceleration  responses 
experienced  by  the  pilot.  The  frequency  and  acceleration  sensitivity  point 
plotted  on  the  MIL-SPEC  requirements  should  not  be  compared  with  lines  of  con¬ 
stant  CAP,  but  with  lines  of  constant  CAP',  which  is  now  not  equal  to  oin^/n/a. 

An  example  will  serve  to  illustrate  this  point. 

Consider  the  A-6  airplane,  including  a  first  order  servo  actuator  lag, 
represented  bv: 
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a  flight  condition  of  .88M  at  20,000  feet  altitude,  the  pitch  rata  transfer 
function  to  elevator  inputs  is  expressed  as: 


ec%) 

SeCs'i 


-  Cs-t  \34p 

S2  +  3.^2  Ms  2  8.£> 


03) 


where 


ojn  =  5.348  rad/sec 
n/a  =87.7  g’s/rad 
0(o)  =  -49.15  rad/sec2 
0ss  =  -2.305  rad/sec 

°zss  =  -64-8^  S 
wn2/n/a  =  .759 

<i(o)/n  =  .759 

zss 


This  condition  shows  perfect  agreement  between  u>n2/n/a  and  CAP. 


Including  the  servoactuator  in  the  system  definition: 


e(s)  __  /-a^.33  \/-H<i.iS(sM.^  \ 

SeCs)  “  U+33.33>/\s2‘^  28.C,  j 


where 


'jJn  =  5.348  rad/sec 
n/a  =  37.7  g’s/rad 
9(o)  =  0 

9gs  =  -2.305  rad/sec 
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n  =  .64.84  g 
zss 

wn2/n/a  =  .759 
0(o)/nZgs  =  0 

Obtaining  the  expression  for  6(t)  and  evaluating  it,  6max  is  found  to  be 

®max  =  -38.05  rad/sec^ 

therefore  0n<j  =  .774 

^max^Zss  =  .587 

2 

For  this  case,  there  is  a  large  discrepancy  (23,eO  between  tun  /n/cx  and  0max/ nzs  • 
However, 

CAP'  -  0„a  =  =  .587 


is  found  to  be  identical  to  the  ratio  of  maximum  pitch  acceleration  to  steady 
state  normal  acceleration.  In  order  to  interpret  this  condition  in  terms  of 
the  MIL-SPEC  requirements,  u>ny  6ncj  must  be  plotted  versus  n/o.  This  point 
can  then  be  correlated  with  boundaries  of  constant  control  anticipation 
parameter. 

Expanding  the  definition  of  the  basic  system  to  include  control  system 
components  such  as  feel  system,  actuator  and/or  feedback  dynamics,  the  general 
pitch  response  to  command  inputs  can  be  defined  as 


eu)  Ke  (S-  v-re^  Kir  (>■>*;') 

Sets')  s (s1-. 2 JP  (s-ePj') 
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where  n  £  m  and  z^,  p.  represent  the  poles  and  zeros  introduced  by  the 
additional  system  components.  For  the  general  case  in  which  n  <  m  (i.e. 

<  #Pj) 
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Following  DiFranco's  development  and  defining  on  attenuation  factor  relating 
the  maximum  pitch  acceleration  with  control  system  dynamics  to  that  of  the 
basic  airframe  yields: 


eco)i? 


•V  '/-r, 
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K  TT 


;^2JS^VS+^SP’  IT  P ; 


li 


substituting  equations  (16)  and  (17)  into  equation  (7)  yields 


CAP  = 


n/0L 


IT  Pj 

K  TT  z 


and  defining  CAP'  =  9n(j  CAP  results  in  - 


CAP' 


o!,,  IT  pj 
*,CL  KTTii 


OT) 


Equation  (19)  is  the  most  general  formulation  of  the  control  anticipation  para¬ 
meter.  It  reduces  to  Bihrle’s  definition  for  cases  in  which  no  higher  order 
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dynamics  are  present  and  to  DiFranco's  formulation  for  cases  in  which  only 
feed  forward  components  are  present.  Examination  of  equation  (19)  again 
demonstrates  that  simply  plotting  the  frequency  of  the  dominant  oscillatory 
root  pair  versus  n/a  for  higher  order  systems  cannot  correlate  directly  with 
the  requirements  of  MIL-F-8785. 

Continuing  with  the  example  aircraft,  the  A-6  longitudinal  control  system 
contains  washed  out  pitch  rate  feedback  of  the  form  .12s/(s+.5).  Including 
this  component  in  the  total  system  definition  results  in  the  following  pitch 
angle  relationship: 


ScC  s') 


-ifese.z  (s-»  _ 


(zo) 


where  un  =  6.821  rad/sec 

n/a  =  37.7  g's/rad 
6 (o)  =  0 

6ss  =  -2.306  rad/sec 
n2sg  =  -64.88  g 
wn^/n/a  =  1.234 
9 (o) / nZss  *  0 

Evaluating  the  pitch  acceleration  time  response,  6max  is  found  to  be  equal  to 
-33.61  rad/sec^. 

Therefore  6ncj  =  .684 

emax/nzSs  =  >518 

In  this  case,  the  discrepancy  between  mn^/n/a  and  &maxJnzss  has  grown  t0  58%. 
Determining  CAP'  by  equation  19,  however,  yields 


cap' 


(33.  S) 


(C84)  =  -5\8 
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which  agrees  with  the  value  obtained  from  the  physical  interpretation  of  the 
control  anticipation  parameter. 

CAP'  BOUNDARY  DETERMINATION 


It  has  been  shown  that  an  expression  can  be  developed  relating  the  short- 
period  frequency  and  acceleration  sensitivity  to  the  control  anticipation  para¬ 
meter  for  any  general  aircraft/control  system  configuration.  It  has  also  been 
shown  that  the  attenuated  control  anticipation  parameter  thus  obtained  is  not 
equal  to  the  lines  of  constant  u^/n/a  defining  the  short  period  boundaries 
of  MIL-F-8785C.  Therefore,  in  order  to  utilize  this  approach  in  evaluating 
aircraft  responses,  the  boundaries  must  be  determined  with  respect  to  the 
attenuated  response  parameters. 

The  MIL-F-8785B  and  C  short  period  requirements  were  established  from 
flight  tests  of  variable  stability  research  vehicles.  Representative  maneuvering 
tasks  were  performed  and  the  short  period  frequency,  damping  and  acceleration 
sensitivity  characteristics  were  evaluated  as  a  function  of  pilot  opinion 
ratings.  A  portion  of  the  data  utilized  in  that  analysis  was  obtained  from 
DiFranco's  report,  reference  (c) ,  with  the  assumption  that  the  feel  systems  effects 
were  negligible.  DiFranco's  data,  first  with  the  assumption  that  feel  system 
effects  are  negligible  and  secondly  with  feel  system  effects  included,  are 
presented  in  Figures  5  and  6,  respectively.  It  should  be  noted  here  that  all 
of  DiFranco's  data  contained  level  1  damping  ratios.  Therefore,  the  pilot 
ratings  can  be  assumed  to  reflect  only  the  frequency/acceleration  sensitivity 
relationship  which,  since  there  was  a  feel  system  included  in  the  test  aircraft, 
should  be  represented  most  accurately  by  CAP'. 

Inclusion  of  feel  system  effects  changes  the  boundaries  primarily  at  the  higher 
values  of  the  control  anticipation  parameter.  In  this  area,  CAP'  is  appreciably 
less  than  CAP,  as  predicted  by  Figure  4.  as  short  period  frequency  increases, 

9nd  ancl  therefore  CAP',  is  decreased.  In  order  to  use  CAP'  as  the  correlating 

parameter,  the  specification  boundaries  must  be  modified  to  reflect  this  vari¬ 
ation.  The  boundaries  suggested  by  DiFranco  are  included  on  Figure  6.  The 
data  are  repeated  in  Figure  7  in  the  format  specified  by  MIL-F-87S5B  and  C. 

From  this  limited  amount  of  data  it  cannot  be  concluded  which  method  is  more 

advantageous.  It  will,  however,  be  shown  that  the  CAP'  formulation  is  prefer¬ 
able  when  considering  equivalent  systems. 

EQUIVALENT  SYSTEM  REPRESENTATION  OF  CAP' 

The  concept  of  equivalent  systems,  reference  (e) ,  seeks  to  identify  an 
equivalent  frequency  and  acceleration  sensitivity  which,  as  closely  as  possible, 
represents  the  higher  order  response.  This  is  accomplished  by  matching  the 
frequency  response  of  the  higher  order  system  (on  a  Bode  plot)  with  that  of 
a  first  order  numerator  over  second  denominator  (short  period  approximate 
transfer  function)  augmented  by  a  time  delay.  The  merits  of  such  a  matching 
procedure  have  been  widely  discussed  and  will  not  be  repeated  here.  What 
is  of  interest,  however,  is  the  correlation  of  the  resulting  equivalent  svstem 
parameters  to  the  higher  order  system  via  the  control  anticipation  parameter. 
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Returning  to  the  first  example  of  the  A-6  airplane,  including  servo 
actuator  lag,  the  following  equivalent  system  is  obtained: 

-.ozfcs 

0 _  -48. Cs-*  e  ^  ^ 

ScCs}  sz  -*■  3.4s  -*■  zQ.W 


where  ion  =  5.302  rad/sec 
e 

n/ae  =  37.7  g's/rad 

9(o)  =  -48.19  rad/sec2 

ogg  =  -2.298  rad/sec 

n_  =  -64.64  g's 
zss 

(a)n2/n/a)e  =  CAPe  =  .746 
0(o)/nzss  =  *746 

There  is  consistency  within  the  j;n2/n/a  and  CAP  parameters  for  the  equiva¬ 
lent  model  (as  would  be  expected  from  the  definition  of  CAP).  However,  the 
control  anticipation  parameter  previously  developed  for  the  higher  order 
system  is  considerably  less  (CAP'  =  .587).  Although  there  is  excellent  agree¬ 
ment  between  the  overall  high  and  low  order  system  pitch  rate  responses, 

(Figure  8),  there  is  a  noticeable  difference  in  the  initial  portion  of  those 
responses.  As  in  the  case  with  and  without  feed  forward  terms,  there  appears 
to  be  an  attenuation  of  the  maximum  pitch  acceleration  response  between  the 
lower  and  higher  order  systems.  Constructing  an  attenuation  factor  of  the 
form 
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I 


multiplying 


this  by  CAPe  yields 


CAPe  0„i 


MQS 


(.2%) 


For  the  A-6  example  - 

CAP^  r  -  .58^ 


which  is  within  .5%  of  CAP'  for  the  higher  order  system.  Therefore,  a  correla¬ 
tion  between  the  higher  order  and  an  equivalent  lower  order  system  can  be 
established  by  plotting  ojny6nd  versus  n/a  for  both  systems. 

The  case  of  the  A-6  with  washed  out  pitch  rate  feedback  yields  the  follow¬ 
ing  equivalent  system: 


✓  ..  n  -.04  2  S 

©IQ  _  -T5.24  (s*  \-3M  Qe _ 

sZ  +  15.23s  47.5^ 


where 

wn 

ue 

=  6.9  rad/sec 

n/ae 

=  37.7  g's/rad 

9  (o) 

=  -75.26  rad/sec“ 

^ss 

=  -2.12  rad/sec 

n  =  -59.65  g' s 
4ss 
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(mn2/n/a)e 

=  1.263 

P(o)/nZss 

=  1.262 

0nde 

=  .447 

cap; 

=  .564 

Again,  CAP  4  is  found  to  provide  very  good  correlation  with  CAP'  (.518)  while 
oi-n'/n/a  is  considerably  different  and  provides  little  information  about  the 
high  order  system  what  remains  to  be  established,  is  whether  or  not  the 
attenuated  control  anticipation  parameters  thus  defined  can  be  correlated  with 
pilot  opinion. 

FLIGHT  TEST  CORRELATION  OF  CAP' 


In  order  to  evaluate  the  acceptability  of  utilizing  the  control  antici¬ 
pation  parameters  for  specifying  the  short  period  dynamics,  two  extensive 
higher  order  system  flight  research  programs  were  analyzed.  Both  flight  pro¬ 
grams  (references  (f)  and  (g))  were  conducted  by  the  CALSPAN  corporation  on 
the  NT-33  aircraft.  The  first  program,  referred  to  as  NEAL-SMITH,  investigated 
the  effects  of  adding  lead/lag  and  pure  lag  control  system  components  to  basic 
aircraft  dynamics  in  maneuvering  flight.  The  second  program,  referred  to  as 
LAHOS,  investigated  the  effects  of  adding  lead/lag,  pure  lag  and  second  order 
control  system  components  to  the  basic  airframe  under  power  approach  flight 
conditions.  Details  of  the  configurations  investigated  are  contained  in  Tables 
I  and  II,  respectively. 

Lower  order  equivalent  systems  were  determined  for  both  programs  using 
the  McAIR  equivalent  system  fitting  routine,  reference  (h).  Equivalent  systems 
for  the  NEAL-SMITH  data  were  determined  by  NAVAIRDEVCEN  while  the  LAHOS 
equivalents  were  obtained  from  reference  (i) .  In  addition,  the  maximum  pitch 
acceleration  was  obtained  from  time  responses  for  each  of  the  configurations. 

The  resulting  high  and  low  order  response  parameters  are  presented  in  Tables 
III  and  IV. 

Before  comparing  the  data  against  the  specif ication  requirements,  it  must 
be  determined  whether  CAP  or  CAP'  is  the  most  appropriate  parameter.  In  order 
to  do  this,  CAP'  was  plotted  against  both  CAPe  and  CAPg  as  shown  in  Figures  9 
and  10.  The  short  period  approximate  control  anticipation  parameter,  CAPe, 
exhibits  large  variations  from  the  higher  order  system's  CAP'.  Including  the 
attenuation  factor  between  the  high  and  low  order  systems  considerably  improves 
the  correlation  as  shown  by  CAP'  versus  CAPg.  Based  on  these  results,  and  the 
fact  that  CAP'  is  directly  related  to  the  aircraft's  time  response,  CAP'  was 
chosen  as  the  correlating  parameter. 

Figures  9  and  10  also  provide  additional  information  concerning  the  effects 
of  freeing  La  in  the  equivalent  system  matching  process.  The  CAPe  -  CAP'  corre¬ 
lation  is  improved  by  freeing  La  at  high  values  of  CAP'  (lead/lag  cases)  while 
at  low  values  of  CAP'  (medium  frequency  lag  cases),  it  is  degraded.  This 
inconsistent  variation  of  CAPg  further  points  out  the  problems  associated  with 
freeing  La  in  the  equivalent  system  matching  process,  as  discussed  in  the  lit¬ 
erature  (e.g.,  references  (e)  and  (i)) . 
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TABLE  I 

NEAL-SMITH  DATA  CONFIGURATION  SUMMARY 
(Table  from  reference  (f)) 


CONTROL 

SYSTEM 

:haracteristics 


( 1)  Numbers/Letters  Indicate  Configurations  Simulated 

(2)  =.75  for  coi  =63,  16;  *5  =.67  for  cl>5  =75 


NOTE: 


man 
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TABLE  II 


LAI  10 S  DATA  CONFIGURATION  SUMMARY 
(Table  from  reference  (g)) 


CONTROL 

SYSTEM 

DYNAMICS 

SHORT  PERIOD  DYNAMICS 
(Nominal) 

Vlnd  -  120  Kt 

r>y/(X  =  4.5  g/rad;  =1.4  sec 

^sp/  sp 

r2 

1.0/. 74 

.  2. 3/. 57  2. 2/. 25 

2.0/1.06 

-  -  -  

3. 9/. 54 

0.4 

0.1 

- 

- 

1  -A 

2-A  | 

0.3 

0<1 

- 

1  -  B 

1 

i 

0.2 

0.1 

- 

1-C 

2-C  j  3-C 

4-C 

E 

0 

- 

- 

1-1 

2-1 

3- 1 (3-0)* 

« 

l 

r—i 

5-1 

0.1 

- 

- 

1 -2 

2-2 

3-2 

0.25 

- 

1-3 

2-3 

3-3 

tom 

5-3 

0.5 

- 

.  :  _  1 

1-4 

2-4 

4-4 

5-4 

1.0 

- 

4-6 

5-5 

5-6 

m 

16/.  7 

1-6 

2-6 

3-6 

■ 

u 

12/.  7 

i 

r 

2-7 

3-7 

4-7 

EEH 

9/.  7 

1-8 

6/.  7 

- 

2-9 

4/. 7 

wmm 

2-10 

4-10 

1 

■ 

c 

16/  .93 

16/ .  38 

1-11 

2-11 

4-11 

5-11 

for  Configuration  3-0  is  2.1/. 12;  for  Configuration  4-0,  2.1/1.23 


CONFIGURATION 

CONTROL  SYSTEM  DYNAMICS 

t^SP  /  **SP 

6-1 

(YF-17  Original) 

('Ss+f)(.  93  s-fl) 

1.9  /.6S 

M 

(.Ss  +  i)(.9Js+r)  (.  o&s *  f) 

1.9/. 65 

(.23*  1)( ■  t  S  +  J)  (  f-  J  5  *■  1  ) 

NOTES:  •  First  number  indicates  base  aircraft  configuration  simulated: 

second  number  or  letter  identifies  control  system  dynamics; 
letters  for  control  system  lead;  numbers  for  lag. 


•  Total  configuration  dynamic  model  includes  feel  system  and 
actuator  dynamics 
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Inspection  of  Tables  III  and  IV  indicates  large  variations  in  pilot  ratings 
for  configurations  possessing  satisfactory,  i.e.,  -  Level  1,  values  of  CAP'. 

For  example,  LAHOS  cases  2-1  through  2-11  all  have  values  of  CAP'  between  .25 
and  1.4  yet  their  pilot  ratings  range  from  2  to  10.  The  equivalent  damping 
ratios  are  all  acceptable  (.35  <  <  1.3)  and  therefore  should  not  contribute 

to  the  variation  in  pilot  rating.  The  parameter  which  has  not  yet  been  con¬ 
sidered  in  this  analysis  is  the  equivalent  time  delay.  This  parameter  is  a 
result  of  the  matching  process.  It  is  used  to  account  for  the  high  frequency 
phase  lags  not  included  in  the  first/second  order  short  period  approximation. 

The  LAHOS  configuration  2  time  delays  vary  from  .067  to  .332  seconds,  spanning 
the  range  of  acceptable/unacceptable  values  specified  in  MIL-F-8785C.  Using 
only  those  configurations  which  exhibit  level  1  damping  ratios,  CAP'  was  plotted 
against  time  delay  for  both  the  NEAL-SMITH  and  LAHOS  data  as  shown  in  Figures 
11  and  12,  respectively.  Simultaneously  applying  the  CAP'  boundaries  shown 
in  Figure  6  and  the  time  delay  requirements  of  MIL-F-8785C  results  in  regions 
of  acceptable/unacceptable  flying  qualities  characteristics  which  agree  reason¬ 
ably  well  with  the  pilot  opinion  data.  The  primary  discrepancy  for  the  MEAL- 
SMITH  data  occurs  in  the  region: 

.25  <  CAP'  <  .6 

.05  <  td  <  .1 

The  level  1  boundaries  could  be  modified  to  exclude  these  data  paints,  however, 
the  boundaries  would  then  be  in  disagreement  with  DiFranco's  dat3  (Figure  6). 
Determining  equivalent  systems  for  DiFranco's  data  results  in  time  delays  of 
.082  -  .085  seconds.  Plotting  that  data  as  a  function  of  CAP'  versus  time 
delay  yields  a  large  number  of  pilot  ratings  of  2  and  3  in  the  region  in  question. 
The  LAHOS  data-  shown  even  better  correlation  than  that  obtained  from  the  MEAL- 
SMITH  data;  the  major  discrepancy  again  occur ing  at  the  lower  values  of  CAP’ 
in  the  level  1  region. 

An  alternate  method  of  analyzing  augmented  aircraft  longitudinal  dynamics 
has  been  proposed  by  Systems  Technology,  Incorporated  (reference  (i)).  It  con¬ 
sists  of  plotting  the  bandwidth  of  the  higher  order  system  (based  on  either 
gain  or  phase  margin)  versus  the  equivalent  time  delay.  STI's  results  for  the 
MEAL-SMITH  and  LAHOS  data  are  presented  in  Figures  13  and  14,  respectively. 

The  same  data  points  which  do  not  correlate  with  the  CAP'/td  boundaries  are 
the  ones  which  determine  the  level  1  bandwidth  boundaries.  Since  there  are 
only  a  limited  number  of  data  points  which  either  (1)  do  not  correlate  with 
existing  boundaries  or  (2)  are  being  used  to  define  new  boundaries,  it  is 
recommended  that  additional  data  be  acquired  and/or  analyzed.  This  analysis 
should  include  a  determination  of  the  relationship  between  CAP'  and  bandwidth. 

The  NEAL-SMITH  data  of  Figure  11  also  indicates  that  the  level  2  boundary 

may  be  extended  to  lower  values  of  CAP'  and  that  it  may  be  possible  to  define 

different  minimum  CAP'  boundaries  for  level  2  and  3  flying  qualities. 

The  data  of  Figures  11  and  12  also  indicate  the  influence  of  damping  ratio 
on  pilot  opinion  for  the  tasks  performed.  In  all  cases,  the  reduction  of 

equivalent  damping  ratio  below  the  minimum  level  1  boundary  (ce  <  .35)  results 

in  at  least  a  one  level  degradation  of  pilot  opinion  (i.e.,  level  1  is  degraded 
to  level  2  and  level  2  to  level  3).  However,  for  equivalent  damping  ratio 
greater  then  the  maximum  level  1  requirement  (;e  >  1.3),  an  improvement  in 
pilot  ratings  is  indicated.  Again,  further  analysis  is  recommended  since  only 
a  minimum  number  of  data  points  are  available  to  support  this  observation. 
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FIGURE  11.  Time  Delay  versus  CAP 
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FIGURE  13.  Bandwidth  Criterion  -  NEAT.-SMITH  Data 
(Reference  (j)  figure) 


NADC-81186-60 


SniiTH  cah of  iu-ra  s'"fl  it- 

<0(jtts:  •  wBu)  ux.UA.Ibfc  i  rtfivo  ibtk-ft'ii  O  1 


BArJiUtaitij  tJflw  (r«Ji//«c.) 

nw^ 


FIGURE  14.  Bandwidth  Criterion  -  LAHOS  Data 
(Reference  (j)  figure) 
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SUMMARY 

A  longitudinal  control  anticipation  parameter  has  been  defined  which 
correlates  high  order  augmented  aircraft  dynamics  with  both  lower  order  equiv¬ 
alent  descriptions  and  pilot  opinion  ratings.  The  approach  has  been  verified 
with  the  data  from  two  extensive  in-flight  flying  qualities  research  programs 
which  span  the  entire  range  of  pilot  opinion  ratings.  Additional  effort  is 
recommended  to  verify  the  minimum  level  1  boundaries  and  to  determine  the 
relationship,  if  any,  between  the  control  anticipation  and  bandwidth  parameters. 
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